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Abstract 

We study the structure and emission of the SS 433 jets in the X-ray emitting region and in the inner 
and hotter portion inside the X-ray emitting region. In order to consider the jet structure from the inner 
to outer regions we develop the hybrid model combining the conical beam and the model beam whose cross 
section grows with the distance more slowly. We find that the jet beams in the inner and hotter portion 
are of two-temperature and emit a large amount of high energy gamma photons. Our analyses suggest 
the thick absorbing envelope to exist in the SS 433 system. Based on our results, we discuss the possible 
acceleration mechanism for the SS 433 jets. 
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1. Introduction 

The galactic source SS 433 is well-known as an unique object exhibiting a jet phenomenon ( for review, see Margon 
1984 and Vermeulen 1989). The moving behavior of Balmer and He I lines has been interpreted successfully in terms 
of the two opposing jets precessing with a period of 164 days. The velocity of matter outflowing in the jet is found 
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to be 0.26c and remarkably stable, where c is the speed of light. The SS 433 system is an eclipsing binary with the 

orbital period of 13 days. The system may consist of an early type star and an accretion disk around a compact 

object. The compact star may be a neutron star or a black hole, although the definite determination has been still 

waited. 

SS 433 has attracted huge attention from the time of the discovery because of its peculiar behavior. A lot of 
observations have been conducted in the radio, optical and X-ray bands. Many theoretical studies have also been 
put forward. In spite of these enormous efforts the most fundamental problems on SS 433 still remain unsolved. 
Those are the nature of the compact star and the acceleration and collimation of the jet beams. 

Recently, the ASCA satellite has detected X-ray lines of various elements, such as Fe, Ni, Mg, Si, S and Ar, from 
the SS 433 jets ( Kotani et al. 1996; Kotani et al. 1997). In the Fe case both the blue-shifted and red-shifted X-ray 
lines, which are emitted respectively from the approaching and receding jets, have been observed at the same time. 
Fitting the analytical jet model to these ASCA data, Kotani et al. (1996, 1997) have examined the properties of the 
X-ray emitting portion of the jet beams. They have derived the important constraints on the jet parameters such 
as the density, temperature and size of the X-ray emitting region and the mass outflow rate. Their analyses have 
also implied the X-ray absorbing gas to exist in the system. 

We study the structure and emission of jet beams in the X-ray emitting region and in the inner and hotter portion 
inside the X-ray emitting region. Based on the results on the jet structure, we discuss the acceleration mechanism 
for the SS 433 jets. 

We present a simplified model for the SS 433 jets in section 2. We derive the analytic solutions for the structure 
and emission of the X-ray jets in section 3. Deriving the analytic solutions, we show that the jet structure is of 
two-temperature in the inner and hotter portion inside the X-ray emitting region in section 4. In the final section we 
discuss the appropriate jet model, the acceleration mechanism and the absorbing envelope based on our analytical 
results on the jet structure and emission. 

2. Model 

The jet may be divided into several characteristic regions depending on the distance from the central engine, the 
inner region, the X-ray emitting region, the optical and radio emission regions, and the outer region. Here, we focus 
on the X-ray emitting region and the hotter region inside of it. 
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No observational evidence is found for the velocity variation along the jets in SS 433. This fact implies that the 

matter entrainment from the ambient medium and other braking processes are unimportant in the region of our 

interest. Hence, we assume the matter velocity constant, v = 0.26c, along the jets. 

The mass outflow rate M in a jet is written as 

M = Spv , (1) 

where r is the distance from the central engine, p the density and S the area of the jet cross section. For simplicity, 
we consider the beam shapes as expressed by 

S = nR 2 oc r n , (2) 

where R is the radius of the jet cross section and n the constant number (Fukue 1987). In the following we examine 
the jet structure and emission especially for the cases of n = 1 and n = 2. The case of n = 2 corresponds to the 
conical beam, while the case of n = 1 to the jet beam whose cross section grows with the distance more slowly than 
that of the conical beam. 

The matter in the jets cools, expanding and emitting radiations. Writing the temperature and radiative loss rate 
respectively as T and A, we have for the energy equation 
3 k dT kT 1 dp A 



2 \irciH dr pmn p dr v 



(3) 



where tuh is the mass of hydrogen, k the Boltzmann constant and p the mean molecular weight of matter. The 
second term on the left hand side of equation (3) and the term on the right hand side represent the adiabatic and 
radiative coolings, respectively. 

We normalize physical quantities with those at a fiducial point r ; x = r/r 0} y = T/T n , z = p/p and / = A/A . 
Here and in the following the subscript is used to denote the physical quantities at r = r . Using these parameters, 
equations (1) and (3) reduce to the dimensionless equations, 

z = x~ n (4) 

and 
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respectively. In equation (5) the free parameter £ is the ratio of the flow time t/o to the radiative cooling time t r0 , 

to = ^ , (6) 

l r Q 

where tf and t r0 arc defined respectively by 

«/. - r i m 

and 

3 kT 

, 2 fj,m H ,„n 

trO-^~ ■ (8) 

Note that the flow time also expresses the adiabatic cooling time due to the plasma expansion. Equation (5) shows 
that the free parameter £ determines the property of solutions and hence the thermal structure of the jet. 

3. One- Temperature Jets 

We consider here the region of the jet beams, where a typical temperature is ~ 10 keV and X-rays are emitted. 
We assume that the jet is optically thin to X-ray photons, satisfying r e — K e pR < 1, where r e is the optical thickness 
and n e is the electron scattering opacity. In this region the ion-electron scattering occurs so frequently that the 
electrons and ions have the same temperature. 

3.1. Analytic Solutions 
When the radiative cooling rate is given by the thermal Bremsstrahlung emission, 

A = A ff = 5.7 x l0 20 pVf ergs g^V 1 , (9) 

the normalized cooling rate is obtained as 

f = zy/y • (10) 

Equation (5) together with equations (4) and (10) can be solved analytically to yield 
^ + § - 1) fen = 1 



(ii) 

h + 1 & (| - M for n = 2 . 
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Note that the solutions (11) tend to the adiabatic cooling formula, y — 1/x 2 / 3 and y — l/x 4 / 3 respectively for n = 1 

and n — 2, as £ approaches to zero. If we intergrate the Bremsstrahlung emission along the jet, we have for the 

X-ray luminosity emitted from the region outside the distance x, 

L(x) = S r poAffo g(x) ergs s _1 , (12) 
where 

3(l + f &)(^-^r)-|eoln(^) forn = l 
9{x) = < (13) 

+ fOT " = 2 • 

In equation (13) x max denotes the outer end of X-ray emitting region. As an x max we choose here the distance 
beyond which the temperature drops below ~ 0.1 keV. 

Adopting the solution for the conical beam (n = 2) at r/r > 1, Kotani et al. (1996) have calculated the line 
X-ray emission from Fe, Ni and other elements and made the detailed comparison with the ASCA observations for 
SS 433. We list the standard set of jet parameters they obtained in Table 1. In Table 1 we also list the jet parameters 
derived from using the solution for the model with n = 1. The solution of n = 1 also reproduces well the observed 
properties, such as the X-ray luminosity ( <~ 10 36 ergs s _1 ) and temperature ( <~ 20 keV). In both solutions the 
fiducial point is chosen to be the base of the X-ray emitting region of the jet. The length of the X-ray emitting region 
is obtained as £x ~ 10 13 cm. The optical thickness across the beam is r e0 ~ 0.1 and hence the jet plasma in the 
X-ray emitting region is transparent, justifying the assumption made in the energy equation. Note that the value 
of the dimensionless parameter £ is ~ 0.1 or less. This fact shows that the temperature structure is determined 
mainly by the adiabatic cooling loss due to the plasma expansion. We notice that the kinetic energy Lk carried by 
the outflowing matter amounts to more than 10 39 ergs s -1 and exceeds the Eddington luminosity for a star of one 
solar mass. 



3.1.1. Jet inside the X-Ray Emitting Region 

The jet radius Rq of the X-ray emitting region and its distance ro from the central engine are at least several 
orders of magnitude larger than the size of the central engine that drives the jet, a neutron star or a black hole. 
Although Kotani et al. (1996, 1997) have truncated the inner portion of the jet, it is quite natural to assume that 
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the jet beam extends more down to the vicinity of the central engine. In the following let us examine the property 

of the jet expected to exist inside the X-ray emitting region. 

We apply the solution (11) to the inner region adjacent to the X-ray emitting region. Following equations (4) 
and (11), the density and temperature of the jet increase with decreasing distance. Hard X and soft gamma rays 
are emitted from the jet. The energy transfer from ions to electrons through the coulomb scattering becomes less 
efficient at higher temperatures, whereas electrons lose energy promptly radiating X and gamma rays and expanding 
adiabatically. The insufficient energy transfer makes the ion and electron temperatures different at T > 10 9 K (see 
the next section for detail). The jet reaches the critical temperature Ti ~ 10 9 K at the distance r\. Here and in the 
following we denote the physical parameters at the innermost region of the one-temperature jet by the subscript 1. 
We list the physical quantities at r = r\ in Table 2. The relativistic Bremsstrahlung emission rate and adiabatic 
cooling rate, which are mentioned in the next section, are used to calculate £i in Table 2. The optical thickness across 
the beam is less than one, although approaches one, and the jet is still optically thin against electron scattering near 
r = r\. Note that the parameter £i in the case of n = 2 (conical case) exceeds, though a little, one. The radiative 
loss as well as the adiabatic loss also contribute significantly to the plasma cooling near the innermost region of 
the one-temperature jet in the conical case. In the n = 1 case the cooling of the jet plasma is still determined by 
the adiabatic loss. The luminosity of X and soft gamma rays emitted from the region at r > r\ is calculated from 
equations (12) and (13) and is listed in Table 2. In the n ~2 case the luminosity of soft gamma rays emitted from 
the region between the distances n and r exceeds the luminosity of X-rays (<~ 10 36 ergs s _1 ) emitted from the 
region outside the distance r approximately by an order of magnitude, while in the n = 1 case the both luminosities 
are comparable. 

4. Two- Temperature Jets 

We have the two-temperature jet inside the distance n, where the ion temperature is higher than the electron 
temperature. Ions transfer energy to electrons through the ion-electron scattering, which in turn lose energy by the 
radiative emission and adiabatic expansion. The distribution of ion temperature along the jet beams can also be 
described by equation (3), taking \x = 1 and substituting the ion energy loss rate A ie for the radiative loss rate A. 
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The ion energy loss rate A ie is given approximately by 

A = A ie = 7T V E , (14) 

2 m H 

where T now expresses the ion temperature and ve is the ion-electron collision frequency. The ion-electron collision 
frequency is written as 

v E = 3.6 x 10 22 -^ s- 1 , (15) 

where T e is the electron temperature. 

We adopt the relativistic Bremsstrahlung emission with the Compton amplification for the radiative energy loss 
of electrons and approximate the emissivity by 

A ff = 4.9 x 10 16 r/pT e ergs g^V 1 , (16) 

where 77 is the Compton amplification factor. Equation (16) can be used approximately even in the optically thick 
case if the photon diffusion time across the jet is shorter than the flow time. Energy loss rate of electrons due to the 
adiabatic expansion can be approximated by 

A ad = ^l . (17) 
m H t f 

The electron temperature is determined from the balance between the energy acquired from ions and the energy lost 
to radiation and adiabatic expansion. When the radiative loss is dominant over the expansion loss (Aff » A a d), the 
electron temperature is expressed as 

T e = 3.8 x ioV 2 /5 T 2/5 K ; ( 18 ) 
whereas when Aff <C A ad , 

T e = 2.0 x 10 5 T 2/5 K . (19) 

In equation (19) the physical parameters for the n = 1 case in Table 1 are used. If we equate the ion and electron 
temperatures in equations (18) and (19), we can derive the critical temperature T\ above which the jet plasma is in 
the two-temperature regime, 

Ti ~ 2 x 10V 2/3 K for A ff > A ad (20) 



8 

and 
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7i ~ 6.9 x 10 8 K for A ff < A ad . (21) 

We choose the position of r ~ r x as the new fiducial point and normalize the physical quantities with those at 
r = 7"i ; x = r/r-i , y = T/Ti , and z = p/pi. The dimensionless equation (5) with the free parameter ^ replaced for 
£o is again derived from the ion energy equation. The parameter £i is now re-defined as the ratio of the flow time 
to the ion energy loss time or the ion-electron scattering time given by 

tu = y E x ■ (22) 

Furthermore, the ion energy loss rate is written in the dimensionless form as 

/ = V /5 • (23) 

We solve equation (5) together with equations (4) and (23) and derive the analytic solution, 
S^f + f ^i(s^r-l) ^ n = l 

(24) 

for n = 2 . 



y 3/5 = 



Using equation (16) and integrating the high energy emission along the jet, we obtain the luminosity of high energy 
photons emitted outside the distance x. The luminosity L(x) is approximately given by 

L(x) - S 1 r 1 p 1 A ffl h(x) ergs s" 1 , (25) 



where 

' T(! + lei) 2/3 ^ for n=l 
h(x) = { (26) 

| Jjs for n = 2 . 

The outer solutions (11) are connected to the inner solutions (24) at r = r\, thereby yielding a break in the 
temperature distribution. As we proceed inwards along the jet, the ion and electron temperatures increase follow- 
ing equations (18), (19) and (24). Moreover, the ion temperature increases faster than the electron temperature, 
rendering the jet of two-temperature. As we move furthermore inwards, we reach the sonic point where the ion 
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sound velocity becomes equal to the flow velocity. Although derived by ignoring the pressure gradient force and 

by assuming the constant flow velocity, the above solution may be used down to the sonic point as long as the 

approximate estimates of the jet parameters are concerned. We list the physical quantities at the sonic point, which 

are denoted with the subscript s, in the Table 3. Here, for simplicity, the compton enhancement factor is set as 

In the conical jet case electrons lose energy mainly by emitting radiations. In Table 3 the radiative cooling is 
adopted for the n = 2 case. The sonic point is located very far from the central engine. The ion temperature is 
approximately an order of magnitude higher than the electron temperature. The thermal equilibrium is not complete 
among ions since the ion-ion scattering time exceeds the flow time. As seen from the values of £i and £ s , the radiative 
cooling as well as the adiabatic cooling contributes significantly to the energy loss of plasma in the two-temperature 
region of the jet. We find that the luminosity of the high energy photons emitted from the region outside the sonic 
point amounts to <~ 2 x 10 39 ergs s _1 , which is approximately comparable to the kinetic luminosity of the outflowing 
matter in the jet beam. 

In the n = 1 case electrons lose energy mainly by adiabatically expanding. In Table 3 the adiabatic cooling 
is adopted for the n = 1 case. The sonic point is located close to the central engine. As in the n = 2 case the 
analytic solution also yields the electron temperature an order of magnitude lower than the ion temperature and 
large luminosity of high energy photons. Note, however, that at the sonic point the radius of the jet cross section 
R s is an order of magnitude larger than the distance r s . Although assumed in the calculation, the uniformity of 
physical parameters over the jet cross section may not be guaranteed when R^> r. Hence, the physical parameters 
at the sonic point, shown in Table 3, may not be reliable. We expect that the analytic solution for the n = 1 
case can be applicable down to the point where R <~ r. We list the physical quantities there in the Table 4. The 
two-temperature solution for the n = 1 case can extend down to the distance r ~ 2.3 x 10 9 cm and cannot reach the 
sonic point keeping R < r. 

Let us consider the hybrid jet model combining the n = 1 and n —2 cases. We adopt the n = 1 case presented 
above to describe the outer part of the jet beyond the distance r ~ 2.3 x 10 9 cm, while for the inner part we adopt 
the n = 2 case (conical jet). The conical jet extends inwards starting from the boundary with the physical quantities 
shown in Table 4. We approximate the electron cooling rate by the radiative loss rate. The density and ion and 
electron temperatures increase with decreasing distance, following equations (4), (24) and (18). We reach the sonic 
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point at the distance r s <~ 9.2 x 10 7 cm and derive physical quantities such as shown in Table 5. Note that the 

sonic point in the hybrid case locates much closer to the central engine compared to the simple conical case shown 

in Table 3. The luminosiy of high energy photons emitted from the region outside the sonic point is also very large 

and amounts to an important fraction of the kinetic luminosity of the outflowing matter. The other properties of 

the two-temperature jet are more or less similar to those of the simple conical case. 

Fukue (1987) already considered the hybrid model for the SS 433 jets. Contrary to our hybrid model, he adopted 

the n = 1 and n = 2 beams respectively for the inner and outer parts of jets. The sonic points lie close to the 

compact objects similarly to our calculation for the n = 1 case (Table 3). The boundary between the n = 1 and 

n = 2 beams may roughly correspond to the base of the X-ray emitting region. If we adopt a standard set of 

parameters for the physical quantities at the base (Kotani ct al 1996), we obtain the radius of the jet cross section 

at the sonic point which exceeds the distance of the sonic point from the compact stars, especially in the neutron 

star case. The uniformity of the physical parameters, mentioned above, may not be guaranteed also in the hybrid 

model by Fukue (1987). 

We note that the electron-positron pairs and appropriate compton amplification in addition to the pressure gradient 
force should be included in the calculation in order to have more accurate estimates on the jet property especially 
around the sonic point. Even if included, the qualitative properties may not deviate significantly from those obtained 
here. 

5. Discussion and Conclusions 

The jet beams in SS 433 may originate near the compact object, or more specifically in the inner region of the 
accretion disk. The X-ray emitting region of the jet is likely to be located at the distance of > 10 11 - 10 12 cm from 
the compact object. The jet beams must be accelerated and collimated to reach the terminal velocity within the 
distance of < 10 11 - 10 12 cm. The mechanisms of acceleration and collimation for the jet beams in SS 433 are poorly 
understood. In general the gas and radiation pressures and magnetic field are considered as a possible driving force 
for the jet acceleration. 

If the shape of the SS 433 jets is simply conical in the X-ray emitting region and also in the inner region adjacent to 
it, we find that the sonic point should be located very far from the central and compact object. The ion temperature 
at the sonic point is very high and close to the maximum attained near the compact object. It is quite difficult to 
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keep the jet matter so hot when the matter outflows from the jet footpoint near the central object to the distant 

sonic point. The jet matter will suffer from the adiabatic cooling and lower the tempetature substantially unless the 

shape of the jet is completely columnar in the region from the footpoint to the sonic point. Let us turn to the jet 

acceletation in this conical case. The free fall velocity at the sonic point due to the gravity of the compact object 

is more than an order of magnitude smaller than the jet velocity. Hence, the acceleration due to the gas pressure 

may not be applicable in this case. The radiative acceleration is also ruled out since the inner and hotter portion 

of the jets is optically thick along the flow direction and hence photons from outside cannot enter into the jet and 

impart momcntun to the matter. The magnetic driving still remains as a possible acceleration mechanism. These 

arguments indicate that the simple conical model for the SS 433 jets is less likely. 

Promising is the hybrid model presented in the previous section. We adopt the model jet with n = 1, where the 
cross section of the jet grows with the distance more slowly than that of the conical jet, in order to describe the outer 
part of the jet including the X-ray emitting region. We use the conical jet to describe the inner and hotter part of the 
jet. This model explains the properties of observed X-rays well. Furthermore, in the hybrid model the sonic point 
can come close to the central object. Hence, the gas pressure as well as the magnetic force may also be considered 
as a possible driving force for the jet acceleration. The radiative acceleration is again ruled out since the inner and 
hotter portion of the jets is optically thick along the flow direction. We note that the inner and hotter portion of the 
jets emit a huge amount of gamma rays. No observational evidence, however, has not been found for an emission 
of such high energy photons. Hence, we conclude that the inner and hotter portions of the jets are hidden by the 
surrounding matter, most probably by the thick envelope and accretion disk. Note that such an envelope, or a dense 
outflowing atmosphere which envelops the system, is also suggested from the substantial fluctuations of the optical 
magnitudes and the infrared spectrum (Vermeulen 1989; Band 1987). 

The collimation of the jet beams is another important issue in the study of the astrophysical jet. The hybrid 
model above indicates that the jet beam may be squeezed and collimated more strongly around the distance of 
<~ 10 9 cm. We expect that at squeezing a part of the kinetic energy of the jet matter may be dissipated into heat 
and reheat the jet matter. A substantial amount of surrounding matter is also required, again, if it is responsible for 
squeezing, from the pressure balance with the ram pressure of the jet matter. If the thick accretion disk channels 
the jet beam and precesses, the precession of the jet would be a natural consequence. The magnetic confinement 
is another possibility for collimation. Our studies show that the magnetic fields, which yield the magnetic pressure 
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comparable to the gas pressure at the sonic point, are of the order of 10 7 G. This field strength may also be used to 

study and set constraints on the properties of the magnetically driven jets. 
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Table 1. Physical quantities at the base of the X-ray emitting region. 



n 


ro 


Ro 


Po 


To 


Co 


M 


L(xo) 


Lk 




(10 10 cm) 


(10 1() cm) 


(Kn n g cm" 3 ) 


(10 8 K) 




(lO-^oy- 1 ) 


(10 36 ergs s" 


: ) (10 39 ergs s" 1 ) 


1 


6.4 


1.2 


4.1 


2.3 


0.05 


2.3 


1.0 


4.4 


2 


52 


2.3 


1.5 


2.3 


0.15 


3.1 


0.97 


5.9 



Table 2. Physical quantities at the boundary inside of which 



the jet plasma is in the two-temperature regime. 



n 


n 


Ri 


Pi 


7i 


Ci L(xi) 




(10 10 cm) 


(10 9 cm) (10" 


1(, g cm- 3 ) 


(10 8 K) 


(10 36 ergs s _1 ) 


1 


1.3 


5.5 


2.0 


6.9 


0.67 1.9 


2 


12 


5.5 


2.6 


20 


1.6 12 



Table 3. Physical quantities at the sonic point. 



n 




Rs 


Ps 


T s 


^ es 




L(x a ) 




(10 6 cm) 


(10 8 cm) (lO"'' 


? g cm" 


- 3 ) (10 n K) 


(10 9 K) 




(10 39 ergs s _1 ) 


1 


4.3 


1.0 


61 


4.4 


9.1 


0.014 


0.084 


2 


1.1 x 10 4 


5.2 


2.9 


4.4 


17 


0.64 


3.8 



Table 4. Physical quantities at the distance where R — r in the n = 1 case. 



r 


J? 


P 


T 


T e 


c 


L(x) 


(10 9 cm) 


(10 9 cm) (10" 


9 g cm- 3 ) 


(10 9 K) 


(10 9 K) 




(10 37 ergss _1 ) 


2.3 


2.3 


1.1 


4.2 


1.4 


0.23 


1.6 
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Table 5. Physical quantities at the sonic point in the hybrid model. 



r 3 


Rs Ps 


T s 


J- es 




L(x a ) 


(10 7 cm) 


(10 7 cm) (l(T 7 gcnr 3 ) 


(10 n K) 


(10 10 K) 




(10 39 ergss _1 ) 


9.2 


9.2 6.9 


4.4 


1.7 


0.13 


1.3 



